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Abstract 



The relativistic Fokkcr-Planck equation, in which the speed of light c appears as a parameter, is considered. It 
is shown that in the limit c — > oo its solutions converge in L^ to solutions of the non- relativistic Fokker-Planck 
equation, uniformly in compact intervals of time. Moreover in the case of spatially homogeneous solutions, 
and provided the temperature of the thermal bath is sufficiently small, exponential trend to equilibrium in 
V^ . L^ is established. The dependence of the rate of convergence on the speed of light is estimated. Finally, it 

O^ ' is proved that exponential convergence to equilibrium for all temperatures holds in a weighted L^ norm. 
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C^ , 1. Introduction 

The Fokker-Planck equation is a widely used model to describe the dynamics of particles undergoing 
diffusion and friction in a surrounding fluid in thermal equilibrium [12l |. For non-relativistic particles with 
mass m > 0, and in suitable physical units, the Fokker-Planck equation is given by 
> , 

2' 9J+p.V,/ = Ap/ + -Vp.(p/), 0=^. (1) 

^T ■ m kl 



ly^ ' Here / — f{t,x,p) > is the one-particle distribution function in phase space; the independent variables 

are the time t > 0, the position a; G M'^ and the momentum p G R^^ of the particles. In the definition 
(•~^ ' of the dimensional constant 9, T is the temperature of the thermal bath and k is Boltzmann's constant. 

The equilibrium state of ^ is given by the Maxwellian distribution, M = exp{—6\p\'^ /{2m,)), up to a 
multiplicative constant that is fixed by the total mass of the system (which is a conserved quantity). 

In this paper we consider a relativistic generalization of ([1]) first introduced in [8| by stochastic calculus 
methods and re-discovered later in [l[ by a different argument (see [9[ for a review on the relativistic theory 
. , of diffusion, as well as the recent paper [11|). In the same physical units used to write ((T)), the relativistic 

j^ ' Fokkcr-Planck equation is given by 



dtf + mC— ■ Vj;/ == ^p^ 



D''dp,f + -p'f 
m 



(2) 



where 



r)'r)J 



D^^ = !!^{S^^ + IP-^), pO ^ ^m^c^ + |p|2 



p" 

with m denoting the rest mass of the particles and c the speed of light. The equilibrium state of ([2]) is given 
by the Jiittner distribution 3 = e^^^P , again up to a multiplicative constant. 
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The purpose of this paper is twofold. First we prove that ([T]) is indeed the correct Newtonian limit of ([3]); 
in particular we show that, as c — ?> oo, solutions of ^ converge in L^ to solutions of ([T]). This provides 
a further justification of (0) as a meaningful relativistic generalization of ([l}. Our second goal is to study 
the trend to equilibrium for solutions of the relativistic Fokker-Planck equation. The latter problem has 
already been considered in [g, where it was shown that solutions of ^ confined in a torus (i.e., x e T^) 
converge exponentially fast in time in the L^ norm to the Jiittner equilibrium, provided the temperature of the 
thermal bath is sufficiently small. In this paper we study the trend to equilibrium for spatially homogeneous 
solutions of (2)). The assumption of spatial homogeneity allows us to derive more accurate estimates on 
the convergence rate. Moreover it will be shown that, at least within the class of spatially homogeneous 
solutions, the small temperature assumption made in [6[ can be (partially) dispensed of. However in order 
to achieve this we have to leave the natural L^ framework and prove exponential convergence in a weighted 
L^ norm. 

The Newtonian limit problem is studied in Section [5] the analysis of the trend to equilibrium is carried 
out in Section [31 

2. Newtonian limit 

The main purpose of this section is to prove the following theorem. 

Theorem 1. Let < /, /c € C^((0, oo) x MP) be solutions of, respectively, eq. ([1]) and eq. ([2]) with initial 
data < f^^,fc^- Assume that f^^{x,p) = 0, for \x\ > R{c), and R{c) growing at most linearly as c ^ oo. 
Assume in addition that 

r^nifin] ■■= f [(1 + br)|V./^'f + (1 + bDlVp/i'f ] dpdx < oo, (3) 

for 7 > 7 and uj > 9. Then ||/^" — /'"|1li ^- =» \\fc{t) — /(OIIl^ ~^ Oj '^•' '^ ~^ °°' uniformly on compact 
intervals of time. 

Throughout the paper we work with smooth solutions of ^ and ([2]) to avoid technical difficulties. 
Moreover the compact support assumption on f™ in the x variable can be removed by adding suitable 
powers of |x| inside the integral (|3]) (which would also allow to treat the general dimension case). We prefer 
to sacrifice the generality of the assumptions for the benefit of a shorter and less technical proof. 

Remark about the notation: In the following, A < B means that there exists a non-decreasing 
function of time (possibly a constant) C(t), independent of c > 1, such that A < C{t)B. Since we are only 
interested in the limiting behavior as c — > oo, the assumption c > 1 is not a restriction. 

Before proving Theorem [U we show that the solution of the relativistic Fokker-Planck equation inherits 
the bound ([3]) on the initial data. 

Lemma 2. // ^ holds, then T^^.y[f] < oo, for all 7,0; > and t > 0. 

Proof. Let u = V xf ■ We compute 

dt f \p\^\u\^dpdx^-2 I \p\'^D'^dp,u-dp,udpdx + l3{3-j) [ |p|"|updpda; 

JRS Jm.0 Jm.B 

+ 1 I dp,{\p\^-h,D'^)\u\''dpdx. 
Since p,D'^ = {mc)-^p>p°, we have dj,j{\p\^-'^p,D'i) < 1 + |p|". Therefore 

dt I \p\^\u\^dpdx< f {l + \p\'^)\u\^dpdx. 



The bound on the integral of |p|'^|wp follows. The estimation for the integral of b|'''|Vp/p is similar, 
although the calculations are more involved. We omit the details. D 



Proof of Theorem]^ The difference Sf= {f — fc) is a smooth sokition of 

dtSf + P ■ VJf - - Vp • {p5f) - ApSf - g, , 



(4) 



where 



gc = Apf,~dp,{D'^dp,fc 



P ■ ^xfc 



We can write down the solution of (U) using the fundamental solution of the operator in the left hand side: 

Sfit,x,p)= J^{t,x,p,y,w)Sf{0,y,w)dwdy+ / F{t - s,x,p,y,w)gc{s,y,w) dw dy ds, 

Jr« Jo Jrs 

where T is given for instance on [10|, Eq. (2.5)]. The proof of the properties of J^ that we use below can also 
be found in 'lO]. In the second term we integrate by parts once in the variable w and obtain 

Sf = J^it,x,p,y,w)Sf{0,y,w)dwdy- / V^J^{t ~ s,x,p,y,w) ■ X{fc)is,y,w) dw dyds 



JRO 



+ T{t- s,x,p,y,w) 



1 



w ■ 'Vyfc{s,y,w)dwdyds, 



(5) 



where X is the vector field X* = du^ — D^^d^ ■ It is easy to show that 

i^(/c)i<c-^knv./.i, i^-i <M!. 

w' c^ 

Using this in ([5]) we obtain 

\\Sfc{t)\\L^ < (/ T{t,x,p,y,w)dpdx]\Sf{0,w)\dwdy 

JR6 \JR6 / 

+ ^ / / I^H^y/cl ( / J^{t - s,x,p,y,w)dpdx] dwdyds 
c Jo JR<5 VJR6 / 

+ ^ / / l'^n^«'/c| / I'^wJ^lit- s,x,p,y,w)dpdx] dwdyds. 
c Jo JR<5 \JR6 / 

Estimating the integrals in the variables {x,p) we obtain 

II'5/W||li<||<5./(0)!Ui + ^ / / \w\'\Vyfc\dwdyds + ^ f -=L 

c Jo JRi3 c jg V^ — 



|wp I Vtu/c| dwdyds. (6) 



By the finite propagation speed property of the relativistic Fokker-Planck equation proved in [l| , and the 
assumption that /"' = for |.t| > R, the solution of 1^ satisfies /c = for |a;| > R + ct. Whence 

w\'^\'^wfc\dwdy < / / \\/^fc\dwdy+ / / |wp|V^/c| dw dy 

J\x\<cJ\w\<l J\x\<cJ\w\>l 

1 /9 / \ -^/-^ 1 /9 

V,„/,|2dwdy) +(/ \w\^-^dw] ( f Iwl^lV.J^fdwdy' 



<r?l-^ 



and for 7 > 7 the integral in the left hand side is Oii:?'^. By exactly the same argument 



/ \wf\^yf,\dwdy<^l'' ( f \Vyf,fdwdy 

JR6 \Jr6 



1/2 



1/2 



f-'^dw 



|u>|>l 



Iwl'^lVyfcl^dwdy 



1/2 



and for w > 9 the integral in the left hand side is 0(c'^/^). Using these estimates in ^ we get 

\\Sf{t)h^<\\Sf{0)h.+O{l/Vc) 
and the theorem follows. 

3 



D 



3. Trend to equilibrium 

In this section wc restrict to spatially homogeneous solutions of ([2]). Moreover for the analysis of the 
trend to equilibrium it is more convenient to rewrite the relativistic Fokker-Planck equation in terms of 
h = f/3. We obtain 



dth = 5„ 



mc /^,^ , pV 



pO y rn^c^ 



\ 10 

;]dr,jh p- \7„h, or equivalently, dth = Ai^'>h + Wh, (7) 

■/ J m ^ 

where the Riemannian metric g and the vector field W are given by 

5., = — fp%,-^V Wh = W^d,.h, W^^--(e + -^)p^ (8) 

7nc \ P J in V ^p cj 

and Ap ■^ denotes the Laplace-Beltrami operator of the metric g. Note that Wi = gijW^ = dpi logu, where 
u denotes the function 

-Rrrfi I 

■s/deFg Y p" 
Let 

diie ^ Z-^-^^P" dp, Z^ [ e-^^P'dp, (10) 

so that d^g is a probability measure. The reason to emphasize the dependence of the measure /i on the 
parameter 9 will become clear soon. In the following we denote by /i a solution of ([7]) normalized to a 
probability density measure: 



u^l^L^.l^e-'^P". (9) 



\\h\\L^dij.a) = / hd^lg = 1. 

This normalization can always be achieved by rescaling the solution. The entropy functional and the entropy 
dissipation functional are defined by 

©[/i] = / hloghdfig, '^[h] = / g{dph,dp\ogh)d^g, 

and the following entropy identity holds: 

j^DMt) = -3[h]it). (11) 

A solution of (O is said to converge to equilibrium in the entropic sense if S)[/i] — >■ = S)[l] as t — > oo, and 
with exponential rate if 'S[h] = 0(e~^*), as t — >■ oo, for some A > 0. A sufficient condition for exponential 
decay of the entropy is the validity of the following logarithmic Sobolev inequality: 



/ hloghdiig <a g {dph, dp log h)diJ.g, for some a > 

JR3 JR3 



(12) 



and for all sufficiently smooth probability densities measure h (not necessarily solutions of ([7])). In fact 
using (|T2]) in (|TT|) we obtain 

4s)[/i] < ^-T)[h] ^T)[h]< cxp{-t/a). 
dt a 

The Ciszar-KuUback inequality, \\h — l||Li(d/j8) < VSD, see [Tj, implies that h converges to equilibrium in 
L^{d^g) with exponential rate {2a)~^, or equivalently, the solution of ([7]) satisfies 

11/(0 -3M||Li(drt<e-*/(2"), (13) 

where 3ai denotes the Jiittner equilibrium with mass M = \\f\\L^idx)- Clearly, (J13p provides the most natural 
notion of convergence to equilibrium for solutions to the relativistic Fokker-Planck equation. 

Thus the question of exponential trend to equilibrium in L^ has been reduced to prove that (|12p holds. 



Theorem 3. The logarithmic Soholev inequality (|12p holds for 6 > Oq = 2~^ ' f^''' ^ constant a given by 
^" 1 y {^0mc^ + i^^Ad-^m'-c^ - 39) , ife>^, 



where 7 (a;) is the rational function 



nx) ~"^-^ ^"^ ...... ex- mc 



Amcx^ 

Proof. The proof is carried out by using the Bakry-Emery curvature bound condition (SJ, Ij] which states 
that (fT^ holds provided the tensor Ric ~ Ric — Vp log u — called the Bakry-Emery- Ricci tensor — satisfies 

Ric > 5^.9- In the definition of Ric, Ric is the Ricci tensor of g, while u is the function Q. In our case the 
Bakry-Emery-Ricci tensor reads 



Now we use 

n(X X\^ I 

mc \ {p^)"^ J p' 



5(^'^)-!lfl^l'-^)^^l^l^ fo-11^^ 



and so 



Ric(X,X) > 



Amc{p^Y 



{20c{py ~ 13(p°)2 + 20to2^V _ ^2^2) 



g{X,X). 



The function on square brackets is J'(po)- It is easy to show that niin{J'(p''),p*' > mc} is strictly positive if 
and only if > Oq. The value of {2a)~^ is obtained by looking for the minimum of T on [mc, oo). D 

The condition 6 > 9o means that the previous result holds only for small temperatures of the thermal 
bath, since 9 ^ T^^. To prove exponential decay of the entropy for all temperatures one needs to find a 
substitute for the Bakry-Emery curvature bound condition used in the proof of Theorem [S] Although there 
are several criteria in the literature for the validity of logarithmic Sobolev inequalities, we were unable to find 
one that applies in our situation. Thus we proceed by a different approach. Since the following argument is 
independent of the dimension, we consider ([7]) with p G R^. Let us consider, instead of the entropy 23 [/i], 
the new functional £,[h] = ||^|1^2(£;„ )■ Computing the time derivative of £,[h — 1] we obtain 



^£[/i-l](i) = -2 / g{dph,dph)dfie. 
dt Jjjjv 



Thus £[/i — 1] decays exponentially, i.e., ft, ^- 1 in L^{diig) exponentially fast, if we show that the following 
Poincare inequality 

/ {h — l)^dfig < A / g{dph,dph) d^ie, for some A > 0, (14) 

holds for all sufficiently smooth probability densities measure h. The validity of the Poincare inequality (jl4[) 
is equivalent to the existence of a spectral gap for the operator in the right hand side of (O, which will now 
be established by applying a criterion due to Wang, see [13|. To adhere with the notation in [l3[, let us 
rewrite ([7]) in the form 

dth = a''dp.dp,h + Vdp,h, t > 0, p e R^, (15) 



where 



- ^,. + 4;V 5^^ , ""^ - '-pP 



yjm^c^ + IpP V TO^c^y \mc^m^c^ -f \p\^ m 



For r > define 

r[Tr{a(p)) + p ■ b{p)] 1 r i \j / \ ■ r a^-'P^Pj 
7(r = sup , C(r) = / 7 s ds, air) = mt r-^. 

Then by [13|, Th.3.1], the spectral gap for the operator in the right hand side of ([T5|) is strictly positive 
provided there exists a function y g C([l, oo)) such that supj>]^ Gy{t) < oo, where 



^^^^)^^/ 



CM /'"gCMMdsd^. 



a(s) 
Theorem 4. T/ie Poincare inequality (|T4| /loMs /or aH > 0. 
Proof. For eq. ([T5|) the function G(i) is given by 






'^^"^ 2/WA r^-Vm2c2+r2 



Let /3 < 6c and pick j/(t) = %_i . After straightforward estimates we obtain 



,N~1 /■* 8r 



Since lini(_i.oo -fl^) = 0, the result by Wang applies and the theorem is proved. D 

Note: While this paper was being written, we have been informed by J. Angst that he was also able to 
prove the Poincare inequality (|14p and therefore the exponential convergence to equilibrium in L^{d^g) for 
solutions of ([7]). The proof by Angst [2| employs a criterion for the existence of a spectral gap to elliptic 
operators established in [5[ . 
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